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Optimal design problems with 
fractional diffusions 


Eduardo V. Teixeira* & Rafayel Teymurazyard 


Abstract 


In this article we study optimization problems ruled by fractional dif¬ 
fusion operators with volume constraints. By means of penalization tech¬ 
niques we prove existence of solutions. We also show that every solution 
is locally of class C 0,a (optimal regularity), and that the free boundary is 
a C 1,J surface, up to a "H n_1 -negligible set. 

1 Introduction 

The goal of this paper is to establish existence and geometric properties of 
nonlocal optimal design problems. This class of problems arises in the study 
of best insulation devices. Motivations also come from semi-conductor theory, 
plasma physics, flame propagations etc. When taking into account eventual 
turbulence or long-range integrations, the model becomes more accurate when 
ruled by nonlocal operators, such as (—A)“. In particular, the model becomes 
sensible to interior changes in the temperature; local versions of the problem 
can only feel changes on the boundary of the body to be insulated. 

Let us recall that the fractional Laplacian (—A)“ is given by 



where PV is the Cauchy principal value and C Uta is a normalization constant. 
The free boundary optimization problem we study here takes the following for¬ 
mulation: given a smooth domain D C R", a nonnegative function ip: D —> R, 
and a positive number ui > 0, minimize the a-energy 



( 1 . 1 ) 


among competing functions u within the functional set 


K(a,w,<p) := {u G H a (R n ) | u = <p in D and C n ({u > 0 }nL> c ) =w} , (1.2) 
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where C n is the n-dimensional Lebesgue measure, and H a (W l ) is the a-fractional 
Sobolev space (see, for example, 0), i.e. the set of functions u for which 


U H« = 


(1 + |£| 2 “)|M(£)l a < oo. 


By usual methods in the calculus of variations it is hard to perform volume 
preserving perturbations as to derive energy estimates. In turn, proving exis¬ 
tence of a minimizer u as well as regularity properties of u and its free boundary 
d{u > 0} n D is in general difficult tasks from the mathematical view point. 

Local versions of the problem have been well studied in the literature, see 
E El HI sa HI ED] among others. A celebrate approach for tackling problems 
involving volume constrains is based on penalization techniques. The idea is 
to introduce an artificial parameter in the energy functional which charges for 
configurations that exceed the volume budget. For fixed values of the penaliza¬ 
tion parameter, the penalized functional can then be analyzed by free boundary 
variational methods. Nonlocal free boundary variational tools were first intro¬ 
duced in [5]. Thus, the starting point of this current paper is to obtain regularity 
results for minimizers of nonlocal free boundary problems with fixed penalized 
parameter. This requires adjustments of existing free boundary methods from 
[6]; some extra work is needed though. Of course all analytic and geometric 
estimates obtained depend upon the penalized parameter, and they blow-up as 
the penalization term goes to infinity. 

The auxiliary penalization problem we consider here takes the following set¬ 
up. Fix an e > 0, we define the e-energy functional 

Je(v):=[ [ ^ (x) dxd y + (£"«“ > 0} n D °)) ■ (1-3) 

j R" Jr" f ■ y i 

where 

fe{s) 

Define the functional set 

K := {uG H a (R n ) \ u = ip in D} , (1.4) 

and then, the penalized problem becomes 

find u £ K such that J E (u) = inf J E {v). (1.5) 

v£K 

As previously mentioned, the minimization problem m is similar to the 
one treated in [B], However, there’s no such thing as a free lunch - the key and 
in general hard issue, though, is to prove that the aimed volume is attained for 
small (but still positive) values of the penalization parameter. This is only pos¬ 
sible by means of a refined control on the rate between volume decreasing versus 
energy increasing, of competing shapes. The appropriate tool for such a control 


J j(s — u>) for s > w, 
( e(s — u) for s < u>. 
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is the so-called Hadamard’s variational formula, whose local (smooth) version 
is known for over one-hundred years, m- Hence, one of the main difficulties we 
handle in this work is to derive a measure-theoretic, a-fractional Hadamard’s 
formula for domain variations. This is accomplished in section U 

Another difficulty in dealing with nonlocal optimal design problems is the 
lack of local information on the unknowns. To bypass this inconvenience, we will 
make use of the extension property discovered in [7]. We consider the upper 
half-space 

R” +1 = {{x,y) el" x R + }, 

and set /3 = 1 — 2a. For u £ C' 2 (R”) we solve the Dirichlet problem 

-div(/Vw) = 0 in R" +1 , (1.6) 

v(x, 0) = u[x). 

A solution to such a problem can be obtained by convolution with the Poisson 
kernel P n ,a{x,y) of the operator div(y /3 V) in R" +1 , see [7]- We have 

td t \ y^ a 

Pn,a{x,y) = Qn, a ^ 2 +y2 )(n+2a)/2 ’ 

where g n>Q , is such that J P n , a {x, 1) dx = 1. From [7] we also know that 

Theorem 1.1. We have (—A) a u(x) = — lim y^v v (x,y). 

s-> o 

Because of the divergence form of the elliptic operator in (11.61) . a Dirichlet 
integral is available. We also observe that if u solves m and d y u = 0 on the 
hyperplane {y = 0}, one can extend it evenly across {y = 0}, and the new 
equation satisfied by u is 

- div(\y\PVv) = 0 in R" +1 , 

v(x,0) = u(x). 

For any open subset H of R ra+1 , we introduce the weighted Hilbert space 

H\/3,n) := {u £ L 2 (H+); | yf/ 2 Vu £ L 2 (H)}, 
where =flnR" +1 . We set 

J e (u,fi):= [ \yf\Vu\ 2 dxdy + f E {C n {{u > 0}nR”nH)), u£H 1 (/3,n). (1.7) 
Jn 

The study of penalized problem (ED is now replaced by the study of local 
minimizers of I E , i.e. functions u that are in F7 1 (/3, B\) and satisfy 

\/B C B±, Vv £ 7F 1 ( / 9,13) with v = u on dB , I E {u,B) < I e {v,B). (1.8) 

The paper is organized as follows: in section [2] we list few analytic and 
geometric properties we will need for the study of the (fixed) penalized functional 
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0). Still in section [5] we show that a minimum of I e is a-Holder continuous, 
which corresponds to the optimal regularity for the free boundary problem. 
We also show that it is non-degenerate. The constants though do depend on the 
penalized parameter e and they blow-up as e —> 0. In section[3]we obtain measure 
estimates on the free boundary. In section 0] we derive a fractional Hadamard’s 
variational formula and in section [5] we prove that for e > 0 small enough - 
but still positive - the volume constraint is verified. This finally provides the 
existence of a minimum for the fractional optimization problem with volume 
constraint. 


2 Preliminaries 

We start off by listing few properties we will need along the article. Initially, 
let us set, for the sake of notations, that a minimizer of (Ill) will be denoted 
by u instead of u e . Here and afterwards by B r (x,y) we denote the ball in R" +1 
centered at (x,y) and of radius r. When x = 0, we simply write B r . We will 
also write B™(x) for the ball in R” centered at x with radius r. We now prove 
existence of minimizers for the £-penalized problem. 

Lemma 2.1. Given a smooth boundary datum, ip(x,y) defined on dB\ nR" +1 , 
problem cm has an absolute minimum u. 

Proof. Take u 0 with £ n ({u 0 > 0}) < w, then / E (w 0 ) < C (uniformly in e), also 
I e > —u). Therefore, there exists a minimizing sequence {rtfcjfcgN- The sequence is 
bounded in H a (Bi) and, because H a is continuously embedded into B 2n l^ n ~ 2a \ 
the sequence {uk}ke n converges (up to a subsequence) to a function u strongly 
in £ 2ra /( Il - 2a )- e j f or an y g > 0, and almost everywhere in R™. Thus, 

£”({;« > 0}) < liminf C n ({uk > 0}) 

k—too 


and 

[ M^lVu) 2 < liminf f \y\ p \S7u k \ 2 ■ 

Jb x k ^°° Jb x 

Hence, u £ H 1 {fi , B\) and since f E is a continuous and increasing function, one 
has 

h(u) < liminf I s (uk) = inf I s {v). 

fc-> oo VEH 1 (^,Bi) 

Therefore u is an absolute minimizer of I e in H 1 (/3,l?i). □ 

For the behavior of a minimizer in its positivity set, we state the following 
proposition and refer the reader for its proof to Proposition 3.1 of [6]. 

Proposition 2.1. Let u be a local minimizer in (EH; and xo G R n be such 
that u(xo,0) > 0. Then 

lim \yf\/u(x 0 ,y) = 0. 
y^o 


4 


Moreover, if u is defined in R n+1 , is positive outside the hyperplane {y = 0} 
and satisfies div(|y|^Vzt) = 0 in its positivity set, together with the estimate 
u(x, y) = 0(\(x, y)\ a ), then (— A) a u(-, 0) = 0 in R n (~l {u > 0}. 

We now turn our attention to optimal Holder estimates for minimizers. In 
the proof of the following theorem we will use the characterization of Holder 
functions, M- given a £ (0,1), if B is a ball in R" +1 , and if there are C > 0 
and p £ (1, n + 1) such that 

Vz G B, Vr <dist(®,0B), [ \S7u\ p < Cr n+l ~ p+pa , (2.1) 

J B r (x) 

then u £ C°’ a (B). 

Theorem 2.1 (Optimal regularity). If u is a local minimizer of I e posed in Bi, 
then u £ C^“(Bi). 

Proof. For every r £ (0,1) and ( Xo,yo ) £ B\, let us consider the harmonic 
replacement of u in B r (x o, yo) (we have chosen r < 1 — |xo|)j he. the solution of 

- divdyfS/hf 0 ^ 0 ) = 0 in B r (x 0 ,yo), K°’ Vo = u on dB r (x 0 , yo)- (2.2) 

From the translation invariance in x , we may assume Xq = 0. We will use the 
notation h r for the solution of m - Note that u is an admissible Dirichlet 
datum (see Theorems 2.2 and 2.3 of 0). Note also, that for all r > 0 one has 
I e (u,B r ) < I e (h r ,B r ). The latter implies 

f \y\ 0 \Vu\ 2 <[ \yf\Vhr\ 2 + Cr n . (2.3) 

JBr J B r 

Note that although the constant C depends on e, but this does not bother us, 
because e will always be fixed (maybe very small but fixed). 

The rest of the proof is the same as of Theorem 1.1 in [6]. We bring it here 
for reader’s convenience. 

Due to the identity / \y\ l3 Vh r ■ V(u — h r ) = 0, we get from (12.31) 

JB r 

I \yf\V(u-h r )\ 2 <Cr n . 

Jb t 
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Therefore, if r < p < 1 , 


\y\P\Vi 


[ \y\P\V(u~h p + h p )\ 2 
J B r 

\ y f\V(u-h. p )\ 2 + f \y\P\X7h p 
Ip J B r 

< Cp n + 2 / \y\P\Vh p \ 2 

JBr 


< 2 


< Cp n + C[- 


< Cp n + C 


,\ n+1+0 

r 

;) 

/ l^iv/g 2 

7 

J B p 

,\ n+l+0 

r 

;) 

/ Iz/l^iv^l 2 . 

7 

)b p 

<5 fe+1 , // = 

= 8 n , then (12.41) 


(2.4) 


/ \yf\Vu\ 2 <Cp k + Cy6 2 ^ f \yf\Vu\ 2 , 

J B s k -(-1 J B s k 

which in turn, by choosing 5 such that q = C5 2( ^ l ~°^ < 1 and using induction, 
gives 


f \ y \^\^ u \ 2 < t —— 

■>B sk 1 - 9 


Therefore, for all r < 1/2 


\y\ 0 \Vu\ 2 < Cr n . 


(2.5) 


/Br- 


Now, if a < 1/2, then /3 > 0, and 


J B |V«|< \y\-^ ' (j g \y\ 0 \Vu\ 2 ^j ' < Cr n+a , 


which is (ED with p = 1, and so u G C' 0,a (-B 1 / 2 ). 
In case of a > 1/2, we get from (12.511 



|Vu| 2 <r~ 0 



\yf\Vu\ 2 


< Cr n ~P =Cr n ~ 1+2/3 , 


which is (12.11) with p = 2. □ 

Next we prove the non-degeneracy of a minimizer. 

Theorem 2.2 (Non-degeneracy). If u is a local minimizer of (11.811 . then there 
exists a constant Co > 0 such that for all x € B(/ 2 (0) (~l {u > 0}, 

u(x, 0) > Codist(x,d{u > 0})“. 
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Proof. Assume that the conclusion of the theorem is not true. It means that 
for every Co > 0 there is a point 0 £ B", 2 (0) fl {u > 0} such that u(z,0) < 
Codist(2, d{u > 0})“. We denote by d the distance of the point (2,0) from 
the free boundary. We also define 6 := u(z, 0). The contradictory assumption 
implies that <5 can be made as small as we wish. 

From the Harnack inequality, m, there exists c > 0 such that u < c5 in 
B d (z, 0). Now if 7 is a smooth nonnegative function such that 

l{x,y) = 0 in B d/2 (z,0), 7(2;,?/) = 2c in B 7d/8 (z,0)\B 3d/4 (z,0), 
we define 

v(x,y) := min(u(x, y), Sj(x, y)). 

Note that v £ B d (z, 0)), and v = u at the boundary of the ball, so it is 

an admissible test function for (11.81) . Therefore 

I E {u, B d (z, 0)) < I e (v, B d (z, 0)). (2.6) 

On the other hand, from the definition of v one has 


/ \yf\Vv\ 2 < / \ y f\Vu\ 2 + o(S), 

JB d (z,0) JB d (z, 0 ) 

and since v = 0 in B d / 2 (z,0), then 

C n {{v > 0}) < C n {{u > 0}) - C n (B d/2 (z, 0) n R"), 
which implies that 

f e (C n ({v>0})) <f E {C n ({u>0})). 

Therefore, 

I e {u,B d (z, 0)) > I e (v,B d (z, 0)), 

which contradicts (12.61) . We remark that the constant Co may depend on e. □ 

Now, arguing as in |6l Proposition 3.3], it is possible to show that u is in 
fact strongly non degenerate. 

Lemma 2.2. If u is a local minimizer of S3 in B i, and (0,0) is a free 
boundary point, then there is C > 0 such that for r £ (0,1/2), 


supu > Cr a . 

B? 

As a consequence of the Holder regularity result and the non-degeneracy, we 
get the positive density result below. 

Theorem 2.3 (Positive density). Let (0,0) be a free boundary point. If u is a 
local minimizer of EH in B 1, then there is a constant Cl > 0 such that for 
every r > 0 


£”({?? = 0 }nB”) > Cir", C n ({u > 0}flB r n ) > Cir n . (2.7) 
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Proof. In fact, from the non-degeneracy we know that there is y £ Bf such that 
u{y) > Cr a > 0. By Holder continuity, u > 0 in Bg r (y) for a small 5 > 0, which 
gives us the second estimate of (1271) . 

To prove the first estimate of m, it is enough to consider the case r = 1. 
Assume the contrary. There is a sequence of minimizers Uk defined in B \, such 
that 0 G d{uk > 0} and 


lim C n ({u k = 0}) = 0. (2.8) 

/c—>-+oo 

Recall that Uk is uniformly Holder continuous. We may assume that the sequence 
converges to rt 0 uniformly. Moreover, one has 

[ \yf\Vu 0 \ 2 dxdy < liminf [ \yf\Vu k \ 2 dx dy. 

Jb x fc->.+oo J Bi 

For every v agreeing with Uk on dB\ one has I e (uk,B i) < I e (v,B i). Together 
with (12.81) this implies for every v £ U 1 ( / 9,I3 1 ) which agrees with u 0 on dB\, 

f e (£ n (B?))+ [ \y\@\S7uo\ 2 dx dy < I e (v,Bi) 

J Bi 

<f £ (C n (B?))+ [ \y\P\Vv\ 2 dxdy. 

Jb x 

Therefore, uo minimizes the Dirichlet integral over the unit ball of M”, and as 
such, satisfies div(|y| d Vuo) = 0 in B\. Since Mo(0) = 0 and uq > 0, then the 
strong maximum principle, see |4], implies that uq = 0, which is a contradiction 
on the non-degeneracy property. Once again we remark, that the constant C\ 
may depend on e. □ 


3 Further properties of solutions 

The ultimate goal of this section is to prove that the free boundaries of local 
minimizers of CLHD have local finite parameter. The results in this section are 
the analogue of the ones from [3]. 

Proposition 3.1. For a local minimizer u in 12, := —(—A ) a u is a non¬ 

negative Radon measure with support in fin d{u > 0}. 

Proof. Once again we recall the extension result from |?J (see also Theorem ll.il 
above). The proof is now the same as the one of Remark 4.2 of [3]. □ 

In the spirit of [3] a representation theorem can be proven. It plays an im¬ 
portant role in the study of the free boundary. 

Theorem 3.1 (Representation theorem). If u is a local minimizer in H, then 
1 . 7t” _1 (/C 0 d{u > 0} D R n ) < oo, for every compact set /C C 12. 
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2. There exists a Borel fucntion q £ such that 


fi{u) = q £ U n 1 [d{u > 0 }, 
that is for any if G Cq°(Q) there holds 

- j \yfVu-S7ip = ( ifq £ d'U 
Jn J\u> o> 


n— 1 


3. There exist constants 0 < c < C < oo depending on n, 17 and e such that 
for B r {x) C 17 and x G d{u > 0} one has c < q £ (x) < C and 

cr n ~ x < %” _1 (B r (x) n d{u > 0} fl R”) < Cr n ~ 1 . 

4- For 1-L n ~ l almost everywhere in d{u > 0} an outward normal v = v{x.q) 
is defined and furthermore 

u(x 0 + x) = q £ (x 0 )(x ■ i'(x 0 )) ( f + o(|a:|), os x ->• 0 . 

5. > 0} fl K ra ) \ (9*{w > 0} fl K”)) = 0, where d* is the reduced 
boundary. 

6. The reduced free boundary d*{u > 0} fl K." is locally a C ' 1 ’ 7 surface. 

Proof. The first three assertions of the theorem follow as those of Theorem 4.5 
of [5]. For the first assertion we also refer the reader to Theorem 1.1 of 113- 
Note that SI fl {u > 0} fl R" has finite perimeter, thus, the reduced boundary is 
defined as well as the measure theoretic normal u(x), for i £ 3*{« > 0 } (see, 
for example uni). The 5th assertion of the theorem is a consequence of the 3rd 
one and properties of minimizers proved above (see nni). C 1,7 regularity of the 
reduced free boundary follows by [T5] , see also m- 


The proof of 4 is similar to the corresponding one in [9] (Theorem 5.5), but 
since in our case we are dealing with the fractional Laplacian, some modifications 
need to be done (similar to [B]). 

For a minimizer u we will denote by S7_(u) C M” the set where it is 0, and 
by f2+(«) C 1" its positivity set. We will also use T(u) C M n to denote the free 
boundary of it, and T*(u) - the reduced free boundary. 

The reduced free boundary is the set of points Xq at which the following 
holds (see P~2]) : given the half ball (£?”)+( xq) := {(£ — xq) ■ v > 0} fl B™{x o), 
one has 

lim £ n ((B?)+(x 0 )An + (u)) = Q 

r™ C n (B?(x o)) 

Note that from the uniform density of 17 ± one has, as r —» 0 at the free boundary 
point Xq 

Br(x o) n r *(u) C {|(x - x 0 ) ■ u(x 0 )\ < o(r)}. (3.1) 
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Indeed, if u(x) = 0 for ( x — Xq) ■ v(xo) > Sr, there is 7 > 0 such that 
C n (B$ r (x) n {u = 0}) > 7 Sr n , which implies 


. t r((B?)+(a:o)Afi + W) 

lim mi-7--—-r- 

C n (B?(x 0 )) 


> 76; 


a contradiction. 

The same argument is valid, if x £ f2_(u) is such that (x — xq) ■ v < —Sr. 

In order to prove 4, we remark, that in fact it follows from the fact that 
blow-up limits at regular points are one-dimensional. To prove it, without loss 
of generality we assume that v(x-o) = e n . Let 

u r (x, y) := \u(x 0 + rx, ry) 

be a blow-up of u, and let uq be a blow-up limit. We need to prove that uo(x, 0) = 
where q is a constant. There exists a coordinate system (x ', x n ) centered 
at 0 such that 

• il + (uo) = R" (this is because of (13. HI ). 

• (—A ) a UQ = 0 in Cl+(u). 

Define v(x) = (x n )“. By optimal regularity and non-degeneracy, there are posi¬ 
tive constants C\, C 2 such that 0 < C\v < uq < C 2 V. By applying the oscillation 
lemma (see El or Theorem 2.5 of 0 ), one has that there exists A € (0,1) such 
that for all small enough r > 0, 


uo . . wo 

OSCB r — < Aosc B 2r —• 
v v 

Since Harnack constants are invariant under the blow-up scaling, the oscillation 
lemma holds at every scale, the solutions being global. Thus, one may apply it 
all the way down from a ball of radius 2 N r (N being arbitrary large) to a ball of 
radius r. Therefore, ^ is a constant. Of course, that constant depends on e and 
also on the blow-up point. Although for our further work this fact is not very 
important, but by following the argument of [3], which proves 2nd assertion of 
the theorem, one can see that the constant is actually q e appearing in 2 . □ 

Next, we see that points, at which the free boundary has a tangent ball, are 
regular points. 

Definition 3.1. We say that a free boundary point xo £ T(u) has a tangent 
ball from outside, if there is a ball B C Q-(u) such that xq £ B 0 T(u). A point 
xo £ T(u) is said to be regular, ifT(u) has a tangent hyperplane at xq- 

The following result is from [ 6 |. The proof can also be concluded from the 
4th assertion of the Theorem 13.11 together with Theorem 14. II below. 
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Theorem 3.2. (The free boundary condition). Let xq G r('u) be a regular point. 
There exists a constant A e (a) such that 


lim 

X—>Xq 


u(x) 

((x - X 0 ) ■ v(x 0 )) a + 


A e (a). 


4 A fractional Hadamard formula 

In this section we proof a fractional Hadamard’s variational formula, which will 
provide a rate control upon volume decreasing versus energy increasing, for 
competing shapes. 

Let u be a local minimizer. We recall the notation T*(u) for its reduced free 
boundary. For given two points Xi, X 2 G r*(u), the idea is to make an inward 
perturbation around Xi, and outward perturbation around xi in such a way, 
that we do not change very much the original volume, and then compare the 
optimal configuration to the perturbated one in terms of the functional J e . We 
proceed as follows. 

Let p : R —> R be a nonnegative function from Cg°[0,1] with J p = 1. For 
any r G (0, dlst L ^ X2 '> ) and y > 0, we define 


Pr(x,y) 


0, y) + 7 rp( l{x ^ 1,y)l 0), 
(x, y) - Trp(^~?^)v(x 2,0), 
(x,y) 


when (x,y) G B r (x i,0) 
when (x,y) G B r (x2 ,0) 
elsewhere. 


If v is any vector in R" +1 , direct computations show that in B r (xi,0) 


DP r (x, y) ■v = v + (-1) 1 


+ 1 / ,( \{x ~ Xj,y)\ \ ((x-Xj,y),v) 

1P V r J \i x — Xi, y)\ 


v{xi, 0), (4.1) 


where (■, •) the the inner product in K . ra+1 and i takes the values 1 and 2. Note, 
that if 7 is small enough, then P r is a diffeomorphism that maps B r (xi , 0) onto 
itself. Indeed, if 7 sup p'(t) < 1, then P r is a local injective diffeomorphism. 
te[o,i] 

Now, if 7 p(t) < 1 — t, for t G [0,1], 


I Pr(x,y) - {xi, 0)| < |(x - Xi,y) | +7 rp 


\(x-Xj,y)\ j < r 


for any (x, y) G B r (xi,0). Finally, note that P r = Id on dB r (xt, 0), therefore P r 
has to be onto. 

For each r > 0 small enough, we consider the r-perturbed configuration by 


v r {P r {x,y)) = u(x,y). 

The idea is to compare our optimal configuration {w > 0} to its perturbation 
{iy > 0} in terms of the penalized problem (11.81) . For any r > 0 small enough 


11 









and i = 1, 2, we consider the blow-up sequence u l r : -Bi(O) —>• R given by 


K{x,y) := — u(xi +rx,ry). 

From the blow-up analysis (see m we know that the set B\ n {u l r > 0} fl R" 
approaches to {(&, 0) £ B i, {(x, 0) -v(xi, 0)) < 0}, as r — > 0. In order to compute 
the change on the volume of the perturbation, we use the Change of Variables 
Theorem to obtain 


C n (B r ( Xi , 0) C {v r > 0} n R") 


r n I 

' JB r (s:i,0)n{D r >0}nl” 


dz 


' -BiD{u r (xi-t-ra;,ry)>0}nIR n 


dx dy 


' B\ D {u * (xi +rx,ry ) >0}flR n 


det (DP r (xi + rx, ry )) dx dy 


(4.2) 


i+(-irv(i(^)i) 


v(xi, 0)\ dx dy, 

\{x,y)\ / 


as r —> 0. Note that there is a constant C{p) such that for any unit vector v £ S 


C(p) = f 

Jb 1 c 


P'i\( x , y )\)( 1(^)1 ’^) dxdy ~ 


in{((x,o),y><o}nK'‘ 

A similar computation shows that 

£"(B r (® j) o)n{u>o}nR n ) r 

rU JBin{<(x,o)i/(xi,o)<o)}nR" 

as r —> 0. From (TOl) . 6D and m we get 


dx, 


(4.3) 


(4.4) 


1 


£ n (B r (xi, o) n {v r > 0} n R") - r (B r (xi , o) n {u > 0} n R") 


as r —> 0. From the Lipschitz continuity of the penalization f E , we obtain 

f e (c n (B r (xi, 0)n{v r >0}nr)) - f e (c n (B r { Xi , o)n{«>o}nr)) 

< -o(r n ). (4.5) 

£ 

Now we should check what happens with the integral part of the functional. 
Initially we observe that 


rU JB r (x U 0 ) 


'Bin{i 4 >o} 


\y\ l3 \V u (x,y)\ 2 dxdy = 
\y\P\Vu l r (x, y)\ 2 dx dy, 


\y\P\^7 u r(x,y)\ 2 dxdy 


(4.6) 
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once r*«) is smooth. Next we apply twice the Change of Variables Theorem 
and take into account that P r maps B r (xi 1 0) diffeomorphically onto itself to 
write 


' B r (xi,0 ) 


lyfl'Vvrix^)] 2 dx dy 


4 r 

r ” JB r (xi, 0 ) 

4 / 

r ” JB r {xi, 0 ) 

/ 

JBin{u*>0} 


|y| / 3 |L»P r (P I T 1 (x,y ))^ 1 • Vm(P,T 1 (x, y)\ 2 dx dy 

\yf\ DP r (z,y)~ 1 ■ S7u(z, y)\ 2 \ det (DP r {z,y))\dzdy (4.7) 

\y\P\DP r (xi + rh, ry ) _1 ■ Vu l r (h, y )\ 2 \det ( DP r ( Xi + rh,ry))\ dh dy. 


From (14.11) . using the fact that for any matrix A with |A| < 1, 


(id + A)- 1 = id + ^-iyA*, 

i =1 

we have 


DP r (xi + rh, ry) ■Vu l r (h,y) 


(4.8) 


= Vu' r (h,y)-{- l ) t+i 7 


On the other hand, 


i + xjm,y)\) 


l(M) I 


((h, y), S7u l r {hi y))u(xi,0) + 0 ( 7 ). 


\det (DP r (xi + rh,ry)) \ = 1 + (~1)* +1 7 P (( h > v)> v i x u °))~ ( 4 - 9 ) 

Combining (14.61) . (14.71) . (14.81) and (14.91) . we obtain 

4 [ \y\ l3 [\S7v r {x,y)\ 2 -\X7u{x,y)\ 2 ]dxdy 

r J B r {xi,0 ) 

f \y\ 0 \\/u l r {h,y)\ 2 p ^ l ' V ^\ {h 1 y),v{x i ,Q))dhdy (4.10) 

jBiniu* > 0 } KM) I 

f y), y))(Vu l r (h , y), v(xi, 0)) d/idy 

JBiflfu’ >0} 

+ 0 ( 7 ). 


= (-l) l+1 7 
+ (— 1 )’ 2 7 1 


l(M) I 


Now we recall that in the proof of 4 of Theorem 13.11 we verified that w* (/i, y) = 
q s (xi)(h n )°[_ + o(r), as r —> 0, i.e. blow-up limits are one dimensional. Therefore, 
from the blow-up analysis, we have, as r —> 0 


l b in{uj.>o} 


\yf\WK \ 2 ^ a 2 qi( Xi ) 




/ B\C\{{z ,Vi) 


\yf\h 


dh dy. 
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Hence, letting r —>■ 0 in (14.101) leads to 


-/ 

r n JB r (xi,o) 

where 


c{p) := lim / 


\y\P [|Vty | 2 - |Vm| 2 ] -+ (-l) i+1 ^a 2 q 2 e ( Xi )c(p) + o( 7 ), (4.11) 

I y\ p P'(\(h,y )|), 


n{w*>o} K\ 1+fl \(h,y)\ 


-{(h,y),v{xi,0)) dhdy 


is a positive constant. To check that indeed, 0 < c(p) < oo is a positive constant, 
we argue as follows: since 

div (p(M)) = p 

so the divergence theorem together with the blow-up analysis provides 

r £%i). 

isin{«*>0} l~l 


(z, u(xi, 0 )) dz —» f p(\z\)dH n 1 (z) = const. 

J Bi n| 


' .Bin{(,z,i/(:ri,0))=0} 

Recalling that the function p is compactly supported in [0,1], we conclude, that 
0 < c{p) < oo. Returning to (14.111) . we can write 


J M^[|Vu r | 2 - |Vu| 2 ] dx dy = r n 'ya 2 c(p) (g 2 (a:i) - q 2 e {x 2 )) + r n o( 7 ). (4.12) 
n 

Combining (14.121) with the minimality property of u, we get 

0 < I e (v r ) - 4 (w) < r n ^a 2 c{p)(q 2 e {x 1 ) - q%(x 2 )) 

+ r n o(j) +-o(r n ), (4-13) 

£ 

which gives after dividing by r n and letting r —> 0 

0 < 'ya 2 c(p) (q 2 {xi) — q 2 (x 2 )) + 0 ( 7 ). 

Now dividing by 7 and letting 7 —*■ 0, and afterwards reversing the places of X\ 
and x 2 , we obtain 

q e {x 1 ) = q e {x 2 ). 

Since 27 and x 2 were arbitrary in T*(u), we actually proved 
Theorem 4.1. On the reduced free boundary, we have q e = A e (a). 

Note also that (14.101) provides the Hadamard’s variational formula: 

[ \y\ P [I Vu r | 2 - |Vu| 2 ] = \ 2 e {a)V + o(V), (4.14) 

J n 

where V is the volume change. 
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5 Recovering the original problem 

Here we shall relate a solution to the penalized problem to a (possible) solution 
of our original problem. The idea is that the function f e will charge a lot for those 
configurations that have a volume bigger than to. We hope that if the charge 
is too big optimal configuration of problem (If .81) will prefer to have volume w 
than pay for the penalization. 

Proposition 5.1. There exist C, c > 0 constants such that 

o < c < c n ({u e >o}nmn)<u + Ce, 

where u e is a solution to Mi- 

Proof. Let f1* be any smooth domain such that its complement contains f l c 
with n R") \ (H c D R")) = w. From the minimality of u e we have 

i e (u e ,n)= [ \yf\Vu e \ 2 + f e (c n {{u e >o}nR n nn))<i e {u„n,) = c, ( 5 . 1 ) 

Jn 

where it* is the a-harmonic function in \ fl c taking ” boundary data” equal 
to p in fl c and 0 on dfl^. Thus 

-{c n {{u e >o}nR"nn)-u)< M£ n {{u. >o}nrnn)<c. 

e 

This proves the estimate from above. In order to prove the estimate from below, 
we first note that since the weight \y\P is in the second Muckenhoupt class A 2 
for ft £ (—1,1), we have a Poincare inequality (see, for example, [ 6 ]), which 
together with (ELD provides 


/ \yf[\Vu e \ 2 + \u E \ 2 ]<C, (5.2) 

Jn 

for some C independent of e. Recalling the fact that u £ takes values p outside 
of the domain, where it is a-harmonic, recalling also the mean value inequality 
and (15.21) . we obtain (by integrating along layers with E := Q fl B$(dn)) 

( [ v) <C(S)C n ({u e >0}nR n nE) [ \yf[\Vu e \ 2 + \u e \ 2 } 

\Jan / Je 

the last integral being bounded uniformly in e. Hence, the estimate from below 
is proved. □ 

Lemma 5.1. There exists C > 0 depending on the domain and p, but indepen¬ 
dent of e, such that X e (a) < C. 

The proof of this lemma is a consequence of Proposition EH isoperimetric 
inequality and 2nd assertion of the Theorem 13.11 (see the proof of the corre¬ 
sponding result in 0, 0>or (23). 
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Lemma 5.2. There exists c > 0 depending on the domain and ip but indepen¬ 
dent of e such that c < A e (a). 

Proof. As in the above mentioned references, the proof is based on the Hopf’s 
Lemma, which is true also for the fractional Laplacian case (see [J or Proposition 
2.7 in 0). 

In fact, let z\ £ fl be such that u e {z{) > 0 for all e > 0. Let S := dist(zi, dfl). 
We then consider the smooth family of domains T t := Bs +t (zi ) fl fl. Let t e be 
the first t such that Ih touches d{u e > 0}. Let xq be that touching point. Define 
if e to be a-harmonic in Tt e \ To, with the following ’’boundary data”: 

ipe 1 3 r 0 = min tp and ip E | r ^ = 0. 

By the maximum principle we have u E > ip E in Tt e \ Tq. From the (generalized) 
Hopf’s Lemma we also know that there exists a constant c > 0 depending only 
on the domain and ip, but independent of e, such that 

ipe(x) > c((x - x 0 ) ■ x 0 )) a . (5.3) 

On the other hand we have the following asymptotic development around Xq 

ife{x) < u E ( x) = A e(a)((x - x 0 ) ■ v(x 0 )) a + o(\x - x 0 |). (5.4) 

Letting x —> xq in (15.411 and taking into account (15.311 , we obtain 

c < A e (a) 


as desired. □ 

Now we are ready to prove the main theorem of this section. 

Theorem 5.1. If £ is small enough, then any solution of (El is a solution of 

( 1 . 1 ). 

Proof. Basically, we just need to show that S := C n ({u e > 0} D R" fl!l) = w, 
for £ small enough. 

Suppose S > uj. In the spirit of the previous section, we consider an inward 
perturbation of the positivity set of u E with the volume change V, such that for 
the new function u E we still have C n ({u E > 0} fl R™ fl fl) > w. Thus 

f E (c n ({u e > o}nr nn)) - f E (c n {{u E > ojnrnsi)) = ~v. (5.5) 

From Theorem 14.11 and Lemma 15. II we have 

/ \yf[\S7h E \ 2 - |Vn e | 2 ] = A 2 e (a)V + o(V) 

Jn 

< C 2 V + o(V). 
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(5.6) 


Using the fact that I e (u e , f2) < I e {u e , O). from (15.51) and (15.61) we get 

0 < C 2 V + o(V) - -V, 

£ 

therefore (by dividing on V and letting V —> 0) it provides us with 

1 

£ - C2’ 

which is a contradiction, when e is small enough. Hence S < to, for small e. If 
S < ui, arguing the same way and using Lemma [5.21 we obtain another lower 
bound for e. Thus, when e is small enough, S = ui. □ 
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